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Preface 


The  purpose  of  this  study  was  to  develop  a  system  to  measure  the 
permittivity  and  permeability  of  Radar  Absorbing  Material  at  millimeter 
wavelengths.  The  Avionics  Laboratory  (AFWAL/AAWP-3)  currently  has  the 
capability  to  make  material  measurements  at  microwave  frequencies,  yet 
modern  radar  applications  may  operate  well  into  the  millimeter-wave 
band. 

It  was  intended  to  test  the  system  at  both  56  GHz  and  94  GHz. 
Although  the  majority  of  experimentation  was  performed  with  the  56  GHz 
equipment  throughout  the  development  of  the  system,  the  56  GHz  oscilla¬ 
tor  failed  before  completion  of  the  project  and  no  final  data  was  taken 
at  56  GHz.  The  absence  of  data  at  56  GHz  makes  it  difficult  to  draw 
conclusions  about  the  performance  of  the  system.  In  addition,  failures 
in  the  receiver  system  forced  modification  of  the  equipment  which  ulti¬ 
mately  reduced  measurement  accuracy. 

I  would  like  to  tnank  the  sponsoring  laboratory,  in  particular 
Edwin  Utt  for  his  contributions  to  the  written  product  and  Dr.  Brian 
Kent  who  conceived  the  project.  I  would  like  to  thank  my  advisor,  Capt 
R.  Jost  for  his  guidance  and  assistance.  Thanks  are  also  due  to  Dr. 

Lair  and  Dr.  Jones  of  the  Air  Force  Institute  of  Technology  for  their 
assistance  with  the  numerical  technique  and  to  the  Fabrication  Shop  for 
the  construction  of  the  sample  holder  assembly.  Most  tf  all,  I  am 
deeply  indebted  to  my  wife  Kathy  for  her  help  and  support.  I  would  also 
like  to  acknowledge  Kathy  and  my  sons  Daniel  and  Scott  for  their 
inspiration,  patience,  and  understanding  throughout  the  course  of  this 
project . 

John  C.  Joseph 
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Abstract 

A  technique  of  measuring  the  permittivity  and  permeability  of 
lossy,  homogeneous,  linear,  isotropic  materials  at  millimeter  wave¬ 
lengths  was  developed.  The  technique  was  tested  by  measurement  of 
representative  materials  at  a  frequency  of  94  GHz. 

For  the  type  of  materials  under  considers  ;ion,  the  permittivity  and 
permeability  are  complex  in  general.  Thus  there  are  four  unknowns  to  be 
determined  and  a  minimum  of  four  measurements  are  required.  The  ap¬ 
proach  taken  was  to  measure  the  magnitude  of  the  transmission  coeffi¬ 
cient  of  a  planar  material  sample  in  free  space  at  multiple  angles  of 
incidence.  The  system  of  equations  obtained  from  any  four  of  the  angles 
is  solved  numerically  for  the  permittivity  and  permeability  of  the 
sample . 

The  equipment  consisted  primarily  of  a  Gunn  phase-locked  milli¬ 
meter-wave  oscillator,  a  pair  of  conical  horn/'lens  antennae,  a  sample 
holder/positioner,  and  a  microwave  receiver.  Harmonic  mixing  was  em¬ 
ployed  to  convert  the  test  waveform  to  the  operating  range  of  the 
receiver . 

For  the  case  of  nonmagnetic  materials,  best  results  were  obtained 
if  the  numerical  solution  was  specialized  to  incorporate  the  assumption 
that  the  permeability  of  the  material  was  equal  to  the  permeability  of 
free  space.  For  a  material  of  arbitrary  permeability,  it  was  necessary 
to  take  data  at  both  perpendicular  and  parallel  polarization  of  the 
incident  wave  to  obtain  satisfactory  results. 


MULTIPLE  ANGLE  OF  INCIDENCE  MEASUREMENT  TECHNIQUE  FOR  THE  PERMITTIVITY 
AND  PERMEABILITY  OF  LOSSY  MATERIALS  AT  MILLIMETER  WAVELENGTHS 


I.  Introduction 


The  effect  that  a  particular  material  will  have  on  electromagnetic 
fields  can  be  predicted  if  two  parameters  of  the  material  are  known,  its 
permittivity  and  its  permeability .  There  are  a  number  of  more 
specialized  and  perhaps  more  familiar  material  parameters;  such  as  index 
of  refraction  in  optics,  intrinsic  impedance  in  microwave  engineering, 
and  power  factor  in  power  engineering;  but  these  can  all  be  expressed  in 
terms  of  £  and  ,  which  are  taken  here  as  the  fundamental  quantities  of 
interest.  Permittivity  and  permeability  are  defined  to  account  for  the 
effects  of  a  variety  of  interactions  between  electromagnetic  forces  and 
the  atoms  of  a  material,  similar  to  the  manner  in  which  modulus  of 
elasticity  is  defined  to  represent  the  elastic  properties  of  a  material 
without  consideration  of  atomic  forces  in  the  study  of  mechanics  of 
materials . 

Permittivity  and  permeability  must  be  obtained  by  experimental 
methods  because  calculation  of  the  material  parameters  is  beyond  current 
theoretical  abilities.  A  technique  for  measuring  £  and  jj[  may  be  based 
on  any  observable  interaction  between  electromagnetic  fields  and  a 
material  sample.  The  choice  of  measurement  approach  depends  on  several 
cons  .derations.  Permittivity  and  permeability  vary  with  frequency  and 
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so  must  be  measured  at  all  frequencies  of  interest.  Although  the  prin¬ 
ciples  are  the  same,  the  nature  of  test  equipm*  '•*  changes  drastically 
with  frequency.  For  example,  bridge  circuits  of  inductors  and  capaci¬ 
tors  may  be  used  up  to  the  megahertz  ranger  while  interferometers  employ¬ 
ing  mirrors  and  light  sources  are  used  at  optical  frequencies.  Another 
consideration  is  the  material  state  (solid,  liquid,  gas)  of  the  samples 
to  be  tested.  It  is  unlikely  that  a  single  sample  holder  will  be  able 
to  accomodate  all  potential  sample  materials,  even  when  restricted  to 
just  one  material  state.  It  should  also  be  mentioned  that  a  significant 
trade-off  exists  between  precision  and  ease/speed  of  measurement. 

The  intended  application  of  this  measurement  facility  is  the  deter¬ 
mination  of  permittivity  and  permeability  of  Radar  Absorbing  Material 
(RAM).  The  present  state  of  manufacturing  of  RAM  is  such  that  two 
samples  taken  from  a  single  sheet  may  have  a  20  to  50  percent  variance 
in  absorber  properties  (1:68).  This  large  variance  has  made  the  goal  of 
extreme  precision  inappropriate  and  has  created  a  need  for  a  system  to 
measure  sheets  of  absorber  quickly  and  without  special  preparation. 

Problem  and  Scope 

The  aim  of  this  thesis  is  the  design,  fabrication,  and  testing  of  a 
system  to  measure  permittivity  and  permeability  at  a  frequency  of  94  GHz 
by  varying  the  angle  of  incidence  of  electromagnetic  energy  illuminating 
a  sample.  Measurement  is  limited  to  solid  materials  in  sheet  form  which 
are  homogeneous,  linear,  and  isotropic.  Multi-layered  sheets  with  dif¬ 
ferent  materials  foi  each  layer  are  not  considered.  No  further  restric¬ 
tions  are  made  to  the  permittivity  and  permeability  of  the  sample  mate 
rial;  measurement  is  not  limited  to  lossless  or  nonmagnetic  materials. 


Assumptions 


In  order  to  determine  permittivity  and  permeability,  one  must 
measure  the  reflection  from  and/or  the  transmission  through  a  material 
sample.  For  this  study,  the  measurement  of  transmission  is  of  interest. 
The  derivation  of  transmission  through  a  sample  assumes  that  an  infinite 
sheet  is  illuminated  by  a  uniform  plane  wave.  The  assumption  of  a  plane 
wave,  which  is  commonly  made  in  electromagnetic  problems,  yet  physically 
unrealizable,  greatly  simplifies  the  analysis.  It  is  prohibitively 
difficult  to  determine  the  form  of  the  actual  field  incident  on  the 
sample  and  oven  more  difficult  to  calculate  the  transmission  of  the 
actual  field  through  the  sample.  A  rigorous  development  of  the  actual 
rather  than  ideal  problem  is  not  warranted  because  the  resulting  com¬ 
plexity  of  measurement  and  increased  computer  requirements  w'"\d  defeat 
the  objective  of  making  measurements  in  a  routine  fashion,  especially  in 
view  of  the  current  manufacturing  tolerances  of  RAM.  The  experimental 
method  of  approximating  the  conditions  assumed  above  will  be  discussed 
in  Chapter  II  of  this  report  and  the  impact  the  assumptions  have  on 
system  design  will  be  seen  in  Chapter  IV. 

General  Approach 

The  transmission  coefficient  was  calculated  for  a  material  of 
arbitrary  £  and  fJi  by  utilizing  a  derivation  (2:87-92)  of  plane  wave 
incidence  on  a  dielectric  sheet.  A  numerical  technique  for  determining 
€  and  fjL  from  measured  values  of  the  transmission  coefficient  at  multiple 
angles  of  incidence  was  devised  and  implemented.  The  equipment  was 
assembled  and  tests  were  made  on  materials  representative  of  pure  die¬ 
lectrics,  nonmagnetic  absorbers,  and  magnetic  absorbers. 


The  information  in  this  thesis  is  presented  in  the  following  order 
Chapter  II  states  the  objectives  and  operating  conditions  of  the  mea¬ 
surement  system.  Candidate  approaches  are  discussed  and  the  rationale 
for  the  approach  selected  is  given.  Chapter  III  explains  the  nature  of 
permittivity  and  permeability  and  derives  the  transmission  coefficient. 
The  derivation  of  several  supporting  relations  is  given  in  Appendix  A. 
The  design  of  the  equipment  set-up  and  the  numerical  solution  to  yield 
6  and ^  is  given  in  Chapter  IV.  Chapter  V  contains  test  data  and  an 
evaluation  of  the  system.  General  conclusions  and  recommendations  for 
further  work  are  given  in  Chapter  VI.  Appendix  B  contains  relations 
needed  for  the  computer  codes. 


II.  Analysis 


Material  measurement  techniques  are  well  established  in  the  direct 
current  to  microwave  frequency  regions  and  at  optical  frequencies. 

These  techniques  are  continually  being  extended  higher  and  lower  in 
frequency.  Due  to  practical  difficulties,  the  millimeter  wave  band 
cannot  be  handled  by  either  microwave  or  optical  techniques  without 
substantial  modification.  The  system  designed  for  this  thesis  more 
closely  resembles  a  microwave  than  optical  system,  but  features  of  both 
approaches  are  incorporated  into  this  design.  This  section  presents  an 
analysis  of  the  design  problem:  it  outlines  several  approaches;  ex¬ 
plains  the  reasoning  for  the  approach  chosen;  and  shows  the  influence  on 
the  design  of  the  type  of  materials  to  be  tested,  the  desired  ease  of 
measurement,  the  desired  accuracy,  and  the  test  equipment  available. 

The  measurement  system  is  designed  to  obtain  £  and  fJi  for  homoge¬ 
neous,  linear,  isotropic  materials.  It  is  intended  to  measure  absorbers 
so  the  assumption  that  the  samples  are  either  lossless  or  nonmagnetic  is 
not  made.  It  will  be  shown  i;i  Chapter  III  that  the  permittivity  and 
permeability  of  such  materials  are  complex  in  general  and  hence  there 
are  four  unknown  quantities  to  be  measured:  th*  real  and  imaginary 
components  of  both  £  and  ^1  • 

The  emphasis  throughout  this  project  has  been  placed  on  simplicity 
of  measurement.  The  degree  of  precision  is  limited  to  the  maximum  that 
can  be  attained  with  existing  equipment  and  without  complicating  the 
measurement  process.  The  test  frequency,  angle  of  incidence,  and  sample 
thickness  are  treated  as  known  quantities  in  the  computations.  It  will 


be  seen  that  the  uncertainty  in  sample  thickness  is  a  large  contributor 
to  the  total  measurement  uncertainty*  Thickness  is  measured  to  a  preci¬ 
sion  of  one  mil  (0.001  inch),  but  variations  in  thickness  of  less  than 
one  mil  can  create  significant  changes  in  the  transmission  coefficient 
at  millimeter-wave  frequencies.  The  wavelength  of  free  space  radiation 
is  approximately  equal  to  126  mils  at  a  frequency  of  94  GHz.  Therefore, 
one  mil  represents  approximately  0.008  wavelengths.  It  is  shown  in 
Appendix  A  that  the  electrical  thickness  (thickness  in  terms  of  the 
wavelength  within  a  material)  depends  on  the  permittivity  and  permeabil¬ 
ity  of  the  material.  The  wavelength  within  a  material  is  necessarily 
smaller  than  the  free  space  wavelength  at  the  same  frequency,  thus  one 
mil  must  represent  more  than  0.008  wavelengths  for  any  material.  A  more 
precise  measurement  of  thickness  would  not  improve  results  because  the 
thickness  of  typical  samples  is  not  uniform  to  within  one  mil.  This  is 
especially  true  for  absorbers  which  have  a  foam  construction. 

The  system  is  also  constrained  by  the  particular  equipment  avail¬ 
able,  a  Hughes  solid-state  continuous  wave  source  of  frequency  94  GHz 
and  a  Scientif ic-Atlanta  phase -amplitude  receiver.  An  approach  is 
sought  that  yields  a  practical  experimental  procedure  as  well  as  an 
analytically  tractable  electromagnetic  problem.  Since  an  exact  solution 
to  plane  wave  incidence  on  a  plane  sheet  is  available,  this  approach  is 
followed  and  an  experimental  technique  that  satisfactorily  approximates 
the  assumed  conditions  is  investigated.  Determination  of  the  four 
unknown  quantities  requires  four  measurements;  the  potential  observables 
are  the  amplitude  and  phase  of  the  reflected  and  transmitted  waves. 

Candidate  systems  for  delivering  an  electromagnetic  wave  to  a 
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material  sample  can  be  classified  as  either  guided  wave  or  free  space 
systems.  Guided  wave  systems  include  hollow  waveguide,  coaxial  line, 
and  various  other  types  of  transmission  lines,  such  as  microstrip  or 
strip-line.  While  these  methods  have  had  success  at  microwave  frequen¬ 
cies,  they  are  unsuitable  in  the  millimeter  wave  band.  All  of  the 
guided  wave  systems  suffer  from  high  transmission  loss  at  millimeter 
wavelengths.  The  signal  amplitudes  are  attenuated  and  measurement  sen¬ 
sitivity  is  reduced,  especially  for  lossy  samples.  Hollow  waveguide  and 
most  transmission  lines  cannot  support  transverse  electromagnetic  (TEM) 
waves,  therefore  they  cannot  propagate  a  plane  wave.  Coaxial  line  can 
support  a  TEM  wave,  but  at  millimeter  wavelengths  the  diameter  of  the 
line  would  have  to  be  very  small  to  avoid  propagation  of  higher  order 
modes.  A  coaxial  sample  would  require  small,  precisely  machined,  toroi¬ 
dal  samples.  Many  materials  of  interest  are  simply  unsuitable  to  be 
prepared  in  this  fashion.  Rectangular  waveguide  would  allow  for  simpler 
sample  shapes,  but  the  small  size  is  still  undesirable  for  measuring  RAM 
because  of  the  inhomogeneities  from  manufacturing. 

The  high  transmission  loss  of  guided  wave  systems  can  be  avoided  by 
use  of  a  free  space  system.  The  free  space  wave  is  a  TEM  wave,  though 
it  has  a  spherical  rather  than  planar  wavefront.  The  usual  method  of 
obtaining  a  good  approximation  to  plane  wave  incidence  on  an  object  is 
to  place  the  object  a  sufficient  distance  from  the  antenna.  The  minimum 
distance  is  chosen  to  set  the  variation  of  phase  of  the  incident  field 
across  the  object  to  a  predetermined  amount.  To  insure  a  phase  varia¬ 
tion  of  no  more  than  7T/8  across  an  object  whose  maximum  dimension  is  D, 
the  minimum  distance  (R)  required  is  given  by  R  **  2D^/^  ,  where  \  is 
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the  free  space  wavelength  (3:936).  Thus  R  is  approximately  12  ft  for  a 
3  in  square  sample  and  48  ft  for  a  6  in  square  sample  at  a  frequency  of 
^4  GHz.  Since  this  condition  would  have  to  be  met  for  the  distance  from 
the  sample  to  the  receive  antenna  as  well  as  the  distance  from  the 
transmit  antenna  to  the  sample,  the  total  distance  required  would  be 
approximately  25  to  100  ft*  In  addition  to  the  obvious  disadvantage  of 
providing  a  large  space  to  house  such  a  system  are  these  disadvantages: 
most  of  the  transmitted  energy  would  not  be  incident  on  the  sample  due 
to  spreading;  the  total  propagation  loss  for  such  distances  would  be 
high;  the  signal  level  available  for  measurement  would  be  exceedingly 
low;  and  it  would  be  necessary  to  construct  an  anechoic  enclosure  to 
shield  against  extraneous  reflections  into  the  signal  path.  If  the 
sample  is  situated  close  to  the  antennae,  the  plane  wave  assumption  is 
violated  but  the  other  problems  associated  with  large  distances  are 
eliminated. 

The  four  measurements  required  may  be  provided  by  placing  the 
sample  normal  to  a  line  joining  two  antennae  directed  at  each  other  and 
measuring  both  the  reflected  and  transmitted  amplitude  and  phase.  This 
arrangement  is  rejected  because  of  multiple  reflections  between  the 
sample  and  antennae.  Also,  the  reflected  wave  may  be  smaller  than 
reflections  in  the  transmit  antenna  and  waveguide  components,  so  the 
measurement  of  reflection  is  unreliable.  An  alternative  to  measuring 
both  reflection  and  transmission  is  to  measure  transmission  alone  under 
two  different  conditions;  such  as  varying  the  frequency,  polarization, 
or  angle  of  incidence  of  the  incident  wave  or  varying  the  thickness  of 
the  sample.  Frequency  must  be  ruled  out  because  the  source  available  is 


not  capable  of  any  variation.  Also,  frequency  is  unsuitable  because  the 
material  parameters  vary  with  frequency.  Using  multiple  thicknesses  is 
undesirable  because  twice  the  number  of  samples  are  required,  the  as¬ 
sumption  that  two  samples  of  the  same  material  have  identical  values  of 
6  and  fji  must  be  made,  and  the  method  relies  too  heavily  on  the  preci¬ 
sion  of  sample  thickness.  Polarization  is  a  good  possibility  if  an  easy 
means  of  changing  it  can  be  arranged.  It  is  undesirable  if  it  requires 
moving  the  antennae  and  waveguide.  The  equipment  should  be  fixed  as 
rigidly  as  possible  since  the  measurements  are  sensitive  to  vibrations. 
Also,  changing  the  equipment  position  to  vary  the  polarization  increases 
the  time  required  for  measurement. 

Varying  the  angle  of  incidence  is  the  best  choice.  The  angle  is 
changed  by  rotating  the  sample  and  keeping  all  other  equipment  fixed.  A 
large  number  of  sets  of  data  can  be  taken  easily  and  quickly.  The 
material  parameters  can  be  calculated  from  each  set  of  data  and  the 
results  averaged  to  eliminate  some  of  the  random  measurement  error.  In 
addition,  the  multiple  reflections  between  the  sample  and  antennae  are 
greatly  reduced  because  the  other  than  normal  incidence  directs  the 
specular  reflection  away  from  the  antennae. 

Since  data  can  be  taken  at  more  than  four  angles  of  incidence,  the 
four  measurements  required  may  be  provided  by  measuring  either  the 
amplitude  or  phase  of  the  transmission  coefficient  at  four  angles  of 
incidence.  It  is  decided  to  measure  amplitude  alone  because  amplitude 
is  simpler  to  measure  and  simpler  to  use  in  the  solution  for  £  and  fjL 
than  is  phase. 

The  measurement  approach  is  to  place  a  planar  sample  between,  but 


close  to,  two  antennae  and  measure  the  amplitude  of  the  wave  transmitted 
through  the  sample  at  multiple  angles  of  incidence.  The  remainder  of 
the  system  is  designed  to  minimize  the  error  from  the  assumption  of 
plane  wave  incidence.  An  infinite  sheet  has  been  assumed  in  the  theo¬ 
retical  calculations  so  that  diffraction  effects  from  the  edges  of  the 
sample  and  scattering  off  the  sample  holder  can  be  neglected.  Experi¬ 
mentally,  this  condition  is  approximated  by  concentrating  the  incident 
energy  to  a  small  portion  of  the  sample  by  the  use  of  narrowbeam  lens 
antennae.  The  lens  antennae  also  serve  to  improve  the  planar  quality  of 
the  incident  phase  front.  It  has  been  shown  that  the  use  of  a  dielec¬ 
tric  lens  can  reduce  the  minimum  range  needed  for  plane  wave  approxima¬ 
tions  (4:252-256).  The  size  of  the  sample  and  the  antenna  to  sample 
distance  are  chosen  to  yield  a  desired  level  of  illumination  on  the 
edges  of  the  sample  according  to  the  particular  antenna  gain  pattern. 

The  incident  wavefront  will  be  spherical  rather  than  purely  planar. 
The  effect  this  has  on  measurement  can  qualitatively  be  compared  to  the 
effect  of  measuring  a  plane  wave  incident  at  a  number  of  slightly  dif¬ 
ferent  angles  of  incidence  and  averaging  the  results  (5:595).  The  error 
introduced  by  the  non-planar  wavefront  can  be  made  negligible  by  choice 
of  the  antenna  to  sample  distances.  This  is  empirically  verified  by 
observing  the  effect  of  varying  these  distances  on  the  measured  value  of 
a  sample  of  plexiglass.  Plexiglass  was  chosen  for  the  verification 
rather  than  a  sample  of  absorber  because  the  material  constants  of 
plexiglass  are  well  known. 


III.  Theory 


The  design  of  a  measurement  system  depends  fundamental! ly  on  the 
nature  of  the  quantities  to  be  measured*  Definitions  of  permittivity 
and  permeability  are  often  given  in  terms  of  a  particular  material 
application,  such  as  use  as  a  dielectric  in  a  capacitor.  The  permit¬ 
tivity  in  this  case  is  derived  from  the  charging  current  drawn  by  the 
capacitor  and  the  loss  current  in  the  material.  For  this  reason,  it  is 
appropriate  to  review  the  theory  of  permittivity  and  permeability  and  to 
state  explicitly  what  is  intended  by  these  terms  in  this  thesis.  This 
section  discusses  permittivity  and  permeability  for  an  arbitrary  mate¬ 
rial  and  then  specializes  to  the  case  of  linear,  isotropic  materials. 

The  discussion  remains  general  regarding  loss  and  the  nature  of  material 
loss  is  described.  A  derivation  of  the  transmission  of  a  plane  wave 
through  a  planar  sheet  is  presented  next.  Subordinate  relationships 
that  are  needed  are  derived  in  Appendix  A. 

The  study  of  electromagnetics  involves  four  vector  field  quanti¬ 
ties:  electric  field  intensity  E,  electric  flux  density  or  electric 
displacement  D,  magnetic  field  intensity  H,  and  magnetic  flux  density  B. 
These  four  quantities  are  often  reduced  to  two  in  practice  by  means  of 
constitutive  relations  between  the  flux  density  vectors  (D  and  B)  and 
the  field  intensity  vectors  (E  and  H)  in  a  particular  medium  or 
material . 

Explanations  of  the  constitutive  relations  are  provided  in  many 
texts  on  electromagnetics.  A  detailed  description  of  the  interactions 
of  electromagnetic  fields  and  materials  down  to  the  atomic  level  is 
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given  by  A.  R.  von  Hippel  (6).  The  information  presented  here  was  taken 
from  von  Hippel,  as  well  as,  Collin  (7)  and  Jones  (8). 


The  constitutive  relations  express  D  and  B  in  terms  of  E  and  H 
respectively  and  the  material  parameters  of  permittivity  and  permeabili¬ 
ty.  These  parameters  are  defined  to  allow  one  to  characterize  the 
effects  of  materials  on  fields  without  consideration  of  the  actual 
mechanisms  acting  on  individual  particles.  That  is,  permittivity  and 
permeability  are  macroscopic  rather  than  microscopic  quantities.  The 
various  interactions  between  particles  and  fields  are  summarized  by  the 
material  properties  of  polarization,  magnetization,  and  conduction. 

When  an  electromagnetic  field  is  applied  t  >  a  material,  forces  are 
exerted  on  the  particles  of  the  material.  The  particles  are  displaced, 
creating  effective  electric  and  magnetic  dipoles  and  possibly  a  current 
within  the  material.  Such  a  current  introduces  loss  since  conduction 
involves  the  transformation  of  electrical  energy  into  heat.  Polariza¬ 
tion  (P)  and  magnetization  (M)  represent  the  effective  electric  and 
magnetic  dipole  moments  per  unit  volume  of  the  material  respectively. 

In  general  £  and  jji  are  complex  tensor  quantities  but  they  reduce  to 
complex  scalars  for  many  materials  and  often  can  be  considered  to  be 
real  scalars  for  practical  purposes.  For  instance,  in  free  space  the 
following  relations  hold: 


Dfs  €o  Efs 


(1) 


Bfs  =  ^  H 


o  “fs 


(2) 


where 
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Farads/meter 


eQ  *  1/ ( 36tt )  x  10'9 

\iQ  =  4v  x  10“  Henrys/meter 

Free  space  is  a  degenerate  case  of  an  isotropic  medium,  a  medium  whose 
properties  are  the  same  in  all  directions.  The  polarization  and  mag¬ 
netization  of  free  space  are  identically  equal  to  zero.  When  an  electro¬ 
magnetic  field  is  applied  to  an  isotropic  medium,  the  resulting  polari¬ 
zation  within  the  medium  is  parallel  to  that  of  the  applied  field.  If 
the  medium  is  linear  and  lossless  as  well  as  isotropic,  the  polarization 
is  directly  proportional  to  the  applied  field.  The  effective  electric 
dipole  moment  represented  by  P  gives  rise  to  a  polarization  current 
density  0P/0t  which  is  a  source  of  magnetic  field.  Maxwell’s  equation 
for  the  curl  of  H  states  that  the  source  of  magnetic  field  in  free  space 
is  the  displacement  current  density,  curl(*H)  *0Dfg^t  .  It  is  con¬ 
venient  to  consider  the  sum  of  the  free  space  displacement  current  densi¬ 
ty  and  the  polarization  current  density  as  a  total  displacement  current 
density.  This  will  allow  an  equation  of  the  form  D  *6  E  to  be  written 
for  any  lossless,  linear,  isotropic  material,  where  P  is  accounted  for 
implicitly  and 

€  =  D/E  (3) 

The  constitutive  relations  for  lossy,  linear,  isotropic  materials 
can  also  be  written  in  the  form  of  Eq  (3)  if  the  vector s*E,  ~H,  "5,  B,  "p, 
and  M  are  replaced  by  their  corresponding  phasors.  When  losses  are 
present,  P  and  M  are  no  longer  in  time  phase  with  E  and  H.  The  use  of 


phasors  allows  one  to  include  the  phase  difference  in  the  constant  of 


proportionality.  Thus  for  any  linear,  isotropic  material 


E  =  ke  I 


M  =  X_  H 

—  m  — 


The  complex  constants  of  proportionality  (k(k  and  km)  lead  to  complex 
values  for  6  and  jj,  . 


Derivation  of  the  Constitutive  Relations 


Combining  the  free  space  displacement  and  polarization  as 


suggested 


e  *  efs  +  e 


For  linear,  isotropic  materials 


2  *  e-0  I  +  Ke  £ 


*  («0  +  V  I 


Permittivity  is  defined  as 


£  =  D  /  E 


=  +  k 

o  e 


^ Vi**  VvVy,*' s V *•’ >  V  V,‘.v  V’.i. \ 


si  ■  »«  **  »  fai  »  ■  ♦  j  *>■  bum  uaai  *•»  km mm *-» * 


The  dielectric  constant  or  relative  permittivity  of  a  material  is 
defined  as 


€  =  €  /  €_ 
j-  '  o 

(ID 

II 

H-* 

+ 

* 

(D 

7  €o 

(12) 

Since  cannot  be  negative,  £r  cannot  be  less  than  one.  Equivalently, 

£  must  be  greater  than  or  equal  to  £Q.  For  lossy  materials,  £  is  com¬ 
plex  and  is  written  as 

€  =  €  *  -  j  €•• 

(13) 

1 

u 

u; 

o 

w 

II 

€"  ) 
r 

(14) 

€  =  €'  -  j  6" 

r  r  J  r 

(15) 

where  €  *  and£"  are  greater  than  zero.  The  minus  sign  is  chosen  in  Eq 

(13)  because  the  imaginary  part  of  £  must  be  negative.  A  positive 

imaginary  part  would  imply  energy  creation  rather  than  energy  loss. 

The  development  for JJi  is  similar  to  that  for  £  with  one  significant 

2 

difference:  the  dimensions  of  P  and  D  are  the  same  (coulombs/meter  ), 
while  the  dimensions  of  M  (amperes/meter)  are  the  same  as  those  of  H, 
not  B.  Thus  an  analog  to  Eq  (6)  cannot  be  written  in  the  form 
B  *  B^s  +  M  .  The  objective  is  to  incorporate  M  into  an  equation  of 
the  form  B  -// H  .  This  can  be  achieved  by  simply  summing  M  and  H. 
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B  =  \i  (  H  +  M  ) 


=  (  H  +  >Sn  H  ) 


=  u  (  1  +  k  )  H 
o  m  — 


Permeability  is  defined  as 


u  =  B  /  H 


*  “o  <  1  +  ) 


Relative  permeability  is  defined  as 


ur  =  u  /  4o 


=  1  +  k 


m 


Similar  to  the  electrical  constants, 


k  >  0 

m  — 


4r  >  1 
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(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 


u  •  -  ^  •- 


-•*  ±  -i  ■  > ■  »  ,Tk  A  .Ti  .~i  ■>  .1 . 1 .1  ■> .i >•>  A  mi  ■"»  .~x  a 


4  =  4'-  j(i" 


4r  =  4  *r  -  J4"r 


Material  losses  can  be  attributed  to  an  actual  finite  conductivity 
as  well  as  to  damping  effects  associated  with  the  various  atomic  inter¬ 
actions.  Allowing  for  a  conduction  current  density  J  and  a  displace¬ 
ment  current  density  Maxwell#s  curl  equation  for  H  is 


curl  (H)  -  jd  D  +  J 


jy€  E  + 


=  ( jd>  +  0 )  E 


=  [jv(€-  -  j€")  +  C j  E 


=  (  j jl r’  +  'ic"  +  ')  E 


where 

cj  *  radian  frequency 
C  *  conductivity 

and  the  substitution  J  m(JE  has  been  made.  Since  the  simple  measure¬ 
ment  of  loss  cannot  indicate  the  nature  of  the  actual  loss  mechanism,  it 
is  convenient  to  represent  all  the  var.  us  losses  by  either  an  effective 
conductivity,  Oeff*  or  an  e^fective  £inary  part  of  the  perrait- 


^  > 


tivity,  0~. 


The  first  approach  yields 


eff 


curl  (H)  =  ( j'x'€  ’  +  aeff)  E 


(33) 


where 


aeff  =  a  +  to 


(34) 


The  permittivity  is  considered  to  be  real  (hence  lossless);  the  conduc¬ 
tion  current  density  “Oeff^  encompasses  the  entire  loss.  The 

second  approach  yields 


curl  (H)  =  jti*  (€•  -  j  €"eff)  E 


(35) 


where 


€"  =  €"  +  a  /  a)  (36) 

The  conduction  :urrent  density  is  considered  to  be  zero;  the  total  loss 
is  attributed  to  a  complex  where 

feff  *  e'  -  j  €"eff  <37> 

Derivation  of  Transmission  Coefficient 

The  approach  taken  is  to  replace  the  configuration  of  a  plane  wave 
obliquely  incident  on  an  infinite  sheet  in  free  space  with  an  equivalent 
transmission  line  circuit  and  then  use  circuit  techniques  to  determine 
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the  transmission  coefficient  (T).  The  derivation  presented  here  follows 
material  given  in  (2:79-92). 

A  cross-sectional  view  of  the  problem  is  shown  in  Figure  1  where: 
d  is  the  thickness  of  the  sample 
0i  is  the  angle  of  incidence 

n  is  the  unit  vector  outwardly  normal  to  the  sample 
k  is  the  wavenumber  vector  of  the  incident  wave 


The  wavenumber  is  a  constant  arising  from  the  solution  of  the  time- 
harmonic  wave  equation.  It  is  equal  to  the  ratio  of  (jJ  and  the  velocity 
of  propagation  of  the  electromagnetic  wave  in  the  medium  (in  this  case, 
free  space).  It  can  be  shown  that  k  is  also  equal  to  27Tdivided  by  the 
wavelength  in  the  medium  (2:21). 


_  A  A  A 

k  =  x  +  y  ky  +  2 

k2  =  Ikl  2 


k2  (38) 

(39) 

(40) 


~  Cl'2  4  € 


(41) 


The  equivalent  circuit  is  shown  in  Figure  2,  where  ZQ  and  Zm  represent 
the  wave  impedance  outside  and  inside  the  sample  sheet  respectively. 
Wave  impedance  in  a  given  direction  is  defined  as  the  ratio  of  the 
electric  and  magnetic  field  components  transverse  to  that  direction. 
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w 


Figure  1.  Geometry  for  Oblique  Incidence 
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Figure  2.  Equivalent  Circuit 

The  wave  impedance  is  a  function  of  and  the  polarization  of 

the  incident  wave. 

For  plane  waves,  E  is  orthogonal  to  k.  Two  cases  of  polarization 
are  considered,  parallel  and  perpendicular.  ine  polarization  is  paral¬ 
lel  when  E  lies  in  the  plane  of  incidence  and  perpendicular  when  is 
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Figure  3.  Normalized  Equivalent  Circuit  with  Waves  Indicated 


Figure  4.  Illustration  of  Waves  at  First  Boundary 

orthogonal  to  the  plane  of  incidence,  where  the  plane  of  incidence  is 
the  plane  containing  n  and  k.  It  will  be  seen  that  T  has  the  same  form 
under  either  polarization,  with  the  expression  for  wave  impedance  corre 


sponding  to  perpendicular  or  parallel  polarization  inserted. 


The  problem  at  hand  is  to  derive  the  ratio  of  transmitted  wave  (a^) 
to  incident  wave  ( a for  the  circuit  of  Figure  3.  The  wave  impedance 
has  been  normalized  to  that  of  free  space  by  defining 

2  *  Zm  '  2o  <«> 

The  electrical  thickness  of  the  sample  is  denoted  by(5*.  The  relation¬ 
ship  between  a^  and  a^  is  due  to  the  effects  of  the  waves  traveling 
through  the  two  boundaries  and  through  the  thickness  (J .  The  behavior 
of  the  fields  at  the  first  boundary  is  illustrated  in  Figure  4,  where 
a^  is  the  complex  amplitude  of  the  wave  indicent  from  the  left 
b2  is  the  complex  amplitude  of  the  wave  incident  from  the  right 
Rmn  is  the  ratio  of  the  amplitudes  of  the  wave  reflected  from  the 

boundary  between  regions  m  and  n  to  the  wave  incident  from  region  m 
Tmn  is  the  ratio  of  amplitudes  of  the  wave  transmitted  into  region  n  to 
the  wave  incident  from  region  m 

From  Figure  4  and  the  definitions  of  R  „  and  T  ^ 

°  mn  mn 


a2  =  Tl2  al  +  R21  b2 


(43) 


ai  =  (a2  -  R21  b2^  /  ti2 


(44) 


& 
'  i* 


bl  =  R12  al  +  T21  b2 


bl  =  R12  ^a2  ~  r21  b2^  ^  Ti2  +  T2i  b 


12  T  1 2 1  2 
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(45) 


(46) 


In  matrix  form 


’al" 

"  1  ~P21  1 

"  a2~ 

1 

"  T12 

1 

-bl- 

-r12  (ti2T21  ”  p12R2l! 

-  b2- 

The  relationship  between  a2  and  a^  is  obtained  by  equating  the 
expressions  for  the  two  rightward  propagating  waves  at  the  right 
boundary : 


a2 


-  a3 


e 


o 


a2  = 


a3 


e 


Similarly 


b2  e 


jt  _ 


b3  e 


b3  e 


-jS 


In  matrix  form 


23 


The  equation  governing  the  behavior  at  the  second  boundary  can  be 


written  by  comparison  to  the  result  for  the  first  boundary. 


V 

1 

1  'R32 

”  a4  ” 

_b3_ 

"  T23 

R23  ^T23T32  ~  R23R32  ^ 

b4 

(53) 


Utilizing  the  definitions  of  reflection  and  transmission  coefficients  in 
terms  of  material  impedances  (2:77) 


Mi 

w 

$sf 

w. 

•Vi 

n*n 


% 
m 

wl| 

«r  •< 


T12T21  ~  R12R21  ^  +  R12^1  "  R12')  *  ^R12^"R12^ 


-  1  -  R122  -  R122 


Similarly,  it  can  be  shown  that 


R23  =  '  R32  =  -  R12 


T23  1  +  R23 


T23T32  R23R32  =  1 


Eqs  (47)  and  (53)  reduce  to 


1 — — 
(N 

m 

i 

cr 

NJ 

1 _ 

m 


Mo 


1  R. 


R23  1 


To  simplify  notation,  denote  and  ^23  by  R.  Substituting  Eqs  (52) 
and  (69)  into  Eq  (68)  and  setting  *  0  (since  there  is  not  a  wave 
incident  from  the  right) 


T12  T23 


1  R  fe 


R  1 


0  e 


1  -R 


L-R  1 


Multiplying  the  matrices  and  substituting 


^12  ^23  ^  +  R)(l  ”  R) 


=  1  -  R' 


yields 


e-^  -R2e--'1  R(e 


R(eJ6-e-J6)  e-J5-R2ej6 


The  overall  transmission  coefficient  T  is  equal  to  a^/a^. 


jS  2  -j6 

e J  -  R  e  J 

a> =  ~ — 34 


T  =  .  i  -  s _ 

e^5  -  R2  e-^6 


(70) 


)1  fa 


J /  VW  VV>,*  V%/Vv%/^  ‘VVvV*  y'vVvVJV  V  V /-‘.V ;V 


^  »  -  » --  ^  *  *•  — 


fc  i  a  *  ^  *  i  *  a a  »  a-a  aa  a  a  sa  autJt  £a£  laama  —- ■  as  --  -  m  **  *«-*  «-■■- -«  »*♦'»  —»  > 


(^i)2 

^Z+l  ’ 


(2=1 ) 2  e-j6 
^z+i' 


(z+1)2  -  (Z-1)2 


(Z+1)2  e^5  -  (Z-1)2  e-^5 


(Z+1)2  e^5  -  (Z-1)2  e~^ 


The  denominator  can  be  simplified  to 


denominator  =  (Z2+2Z+1  )e-^  -  (Z2-2Z+1  )e~-^ 


=  (Z2  +  l )  (e^-e"^' )  +  2Z(e^+e”-^  ) 


=  (Z  +l)[2j  sin(5)j  +  2Z[2  cos(6)] 


4Z  cos ( 6 )  +  2j(Z  +1)  sin(6 ) 


From  Eqs  (78)  and  (82) 


4  Z  cos ( & )  +  2j(Z  +1)  sin ( 5 ) 


cos(5)  +  2 2l£  1 12  sin(S) 

4  Z 


(  34) 


w$m 


✓ 


so 


cos(6)  +  2  (Z+2.)  sin(5) 
2  Z 


(85) 


T  =  [  cos  ( 5  )  +  (j)(Z+i)  sin  ( 6  )  l"*1 

2  z 


(86) 


For  perpendicular  polarization 


Ur  COS  (?i) 

[nr  -  sin2(9i)]^ 


For  parallel  polarization 


[u  9  -  sin2  (  0  .  )  ~Y2 

*-  r  r _  l  _ 

€  cos  ( 9 . ) 

r  l 


(87) 


(88) 


The  electrical  thickness  of  the  sheet  is  the  same  for  both 


6  =  2tt(~)  [|i  €  -  sin2(6.)]'s  (89) 

Ao 

The  derivations  of  Z  and  (5  appear  in  Appendix  A. 

By  utilizing  Eqs  (86),  (87),  (88),  and  (89),  it  is  possible  to 
determine  6  r  and fXx  from  measurement  of  T.  The  method  of  solving  these 
equations  is  shown  in  Chapter  IV. 
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IV.  Design 
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The  basic  approach  taken  to  the  design  of  the  measurement  system 
was  presented  in  Chapter  II*  Due  to  the  limited  power  available  from 
the  millimeter  wave  source  and  the  relatively  high  insertion  loss  as¬ 
sociated  with  millimeter  waveguide  components,  an  approach  was  pursued 
that  led  to  minimum  usage  of  millimeter  wave  equipment.  Measurement  of 
absorber  is  restricted  by  the  sensitivity  of  the  receiver  and  the  power 
level  at  the  input  to  the  receiver.  To  achieve  the  maximum  dynamic 
range  possible,  the  waveguide  loss  must  be  minimized  and  the  test  fre¬ 
quency  must  be  converted  to  the  operating  range  of  the  receiver  as 
efficiently  as  possible.  It  will  be  seen  that  the  simplicity  of  meas¬ 
urement  results  in  complexity  of  the  computer  inversion  of  the  measured 
data  to  yield  6  and  f±  .  This  section  provides  the  details  of  the  design 
and  points  out  the  influence  of  several  practical  considerations  on  the 
design.  The  equipment  set-up  is  discussed  first  followed  by  the  compu¬ 
ter  implementation  of  the  numerical  solution. 

Equipment  Design 

A  simplified  block  diagram  of  the  test  equipment  appears  in  Figure 
5. a.  The  primary  components  are  the  millimeter  wave  oscillator,  the 
receiver  subsystem,  and  the  antenna/sample  holder  assembly.  A  photo¬ 
graph  of  the  antenna  and  sample  holder  is  given  in  Figure  5.b.  The 
assembly  provides  a  rigid  support  for  the  antennae  and  the  sample  and 
allows  accurate  positioning  of  the  sample.  The  assembly  is  constructed 
almost  entirely  of  wood  to  minimize  unwanted  reflections.  Plastic  bolts 
are  used  to  secure  the  sample  to  the  wooden  frame.  The  antenna  supports 
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Figure  5.a  Simplified  Block  Diagram 


Figure  5.b  Antenna/Sample  Holder  Assembly 

are  metal  because  maximum  stability  of  tae  antennae  is  needed  and  refle¬ 
ctions  from  the  antenna  supports  are  not  noticeable. 

Ideally,  it  is  desired  to  transmit  a  plane  wave  onto  a  material 
sample  without  illuminating  external  objects  or  the  sample  holder. 

Since  this  cannot  be  achieved  in  practice,  a  compromise  must  be  made 
that  involves:  the  antenna  characteristics,  the  level  of  illumination 
of  the  sample  holder,  and  the  choice  of  sample  size  (A),  maximum  angle 
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of  incidence  (umax)>  and  antenna  to  sample  distance  (R). 

The  necessity  of  confining  virtually  all  of  the  energy  radiated  by 
the  transmit  antenna  to  a  portion  of  the  sample  requires  the  use  of  a 
narrow  beamwidth  antenna.  The  antenna’s  amplitude  pattern  is  therefore 
strongly  non-planar,  yet  a  nearly  planar  phase  pattern  can  still  be 
achieved  by  fitting  the  antenna  aperture  with  a  properly  shaped  dielec¬ 
tric  lens.  The  effect  of  the  lens  is  to  convert  the  phase  front  from 
spherical  to  planar  as  illustrated  schematically  in  Figure  6.  The 
transmit  and  receive  antennae  are  the  same  model:  a  conical  horn  with  a 
plano-convex  lens.  The  amplitude  pattern  of  the  antenna  is  given  in 
Figure  7. 

The  following  considerations  should  be  kept  in  mind  during  the 
derivation  of  the  relation  between  R,  A,  $max»  and  the  antenna  beam- 
width^  : 

1.  the  larger  R  is,  the  more  planar  is  the  incident  wave  across 
the  sample 

2.  the  choice  of  ^ depends  on  the  maximum  allowable  illumination 
level  of  the  sample  holder 

3.  the  greater  A  is  the  wider  the  range  of  incidence  angles  that  can 
be  used.  It  will  be  seen  that  more  varied  values  of^  are 
desirable  for  the  numerical  solution. 

4.  There  is  a  practical  limit  to  the  sample  size  A. 

As  is  increased  from  normal  incidence  *  0°),  the  antenna  beam- 
widths  subtend  a  larger  are?  on  the  surface  of  the  sample.  This  effect 
is  more  pronounced  in  the  horizontal  plane  than  the  vertical,  since  the 
sample  is  rotated  about  a  vertical  line.  The  antenna  beamwidth.  E-plane 
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Figure  8.  Relationship  between  R,  ,  Omax 


,  and  w. 


the  sample  is  square  and  w  is  sufficiently  large,  the  height  of  the 
sample  is  also  sufficiently  large.  Therefore,  the  derivation  of  the 
relation  governing  in  the  horizontal  and  vertical  planes,  R,  A,  and 
Qmax  can  ke  replaced  with  the  two  dimensional  problem  involving  R,  , 
0roax>  and  w*  Figure  8  shows  the  geometry  of  the  relation  between  R, 

7>8 max>  and  w- 

From  Figure  8  and  the  law  of  sines 
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The  maximum  R  for  a  given  ^(beamwidth  corresponding  to  maximum  allowa¬ 
ble  Illumination  level  on  sample  holder),  w,  and  9max  is  expressed 
by  Eq  (96).  If  R  Is  set,  the  maximum  angle  of  Incidence  for  choice  of  'Y 
and  w  Is 
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To  summarize: 

1.  An  increase  in  R  requires  a  decrease  in  $max  and/or  an  increase  in 

w. 

2.  An  increase  in  0max  requires  a  decrease  in  R  and/or  an 
increase  in  w. 

An  increase  in  w  requires  an  increase  in  R  or  the  planar 
quality  of  the  incident  wave  is  degraded. 

4.  Decreasing  allows  for  larger  $max  ,  larger  R,  and  smaller 
samples . 

The  sample  holder  was  built  to  hold  samples  of  size  up  to  12  in  by 
12  in  and  thickness  up  to  1  1/2  in.  For  a  maximum  sample  holder  illumi¬ 
nation  of  30  db  below  the  main  beam  level,  the  E  and  H-plane  beamwidths 
are  set  to  10°  and  7°  respectively.  Pairs  of  valuer  of  R„„v  and 
are  given  for  both  perpendicular  and  parallel  polarization  in  Table  I. 

In  Chapter  V  the  effect  of  varying  R  on  the  measured  data  will  be  shown. 

A  block  diagram  of  the  entire  equipment  configuration  is  given  in 
Figure  9.  Equipment  failures  in  the  receiving  subsystem  forced  the 
modification  of  the  original  set-up  (Figure  9.a)  to  the  final  equipment 
set-up  (Figure  9.b).  The  primary  component  of  the  system,  the  antenna 
support/sample  holder  assembly,  has  been  described  in  detail.  The 
remainder  of  the  equipment  serves  to  generate  and  receive  the  millimeter 
wave  test  signal. 

The  essential  requirements  cf  the  source  are  frequency  and  power 
stability.  The  source  used  was  a  Gunn  phase-locked  oscillator.  The 
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Rmax  versus  #max 


■Ux  <ln> 

^max 

Perpendicular 

Parallel 

(7=  7°) 

(7=  io°) 

50° 

26.8 

17.3 

55° 

23.1 

14.6 

60° 

19.2 

11.8 

65° 

15.2 

8.9 

70° 

11.1 

6.0 

oscillator  is  locked  to  a  low  frequency  stable  reference  signal  from  a 
self-contained  crystal  oscillator. 

The  receiver  subsystem  employs  harmonic  mixing  to  lower  the  test 
frequency  to  the  normal  operating  range  of  the  receiver.  In  the  origi¬ 
nal  arrangement,  the  eighth  harmonic  of  an  11.45  GHz  local  oscillator 
(LO)  was  mixed  with  the  94  GHz  test  frequency  to  produce  an  intermediate 
frequency  (IF;  of  2.4  GHz.  This  IF  was  filtered  and  amplified  prior  to 
proceeding  to  the  1783  receiver.  In  the  modified  arrangement,  the  test 
signal  had  to  be  mixed  with  a  harmonic  of  the  receiver’s  local  oscilla¬ 
tor.  To  maximize  the  received  power  and  therefore  the  available  dynamic 
range,  it  was  necessary  to  mix  the  94  GHz  signal  with  a  harmonic  as  low 
in  number  as  possible.  The  1783  receiver  was  replaced  with  the  less 
sophisticated  1752  receiver  because  the  1752  had  an  LO  of  1  to  4  GHz, 
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which  allowed  mixing  to  be  accomplished  with  the  twenty-forth  harmonic 
of  a  3.917  GHz  LO  frequency.  The  LO  of  the  1783  was  limited  to  2  GHz, 
so  the  94  GHz  test  signal  had  to  be  mixed  with  the  forty-eighth  harmonic 
of  1.938  GHz.  Although  the  1783  was  able  to  lock  to  the  test  signal, 
the  received  power  levels  were  unsuitably  low.  The  attenuators  were  not 
included  in  the  modified  set-up  because  the  power  levels  were  too  low  to 
saturate  the  mixers  or  receiver.  The  isolators  were  used  to  minimize 
reflections  in  the  waveguide  components.  The  1752  receiver  required  the 
use  of  a  voltmeter  to  obtain  amplitude  data.  Measurement  accuracy 
suffered  from  this  because  as  the  decibel  values  of  power  approached  -20 
to  -30  dB  for  absorber  samples,  the  linear  values  of  voltage  decreased 
by  a  factor  of  10  to  32. 

(T2)dB  ”  10  l°g1o^power  rati°) 

*  20  log^Cvoltage  ratio) 
voltage  ratio  *  log^”1  (T2dB/20) 


Numerical  Solution 

The  transmission  coefficient  (T)  can  be  expressed  as  a  function  of 

and^/M  with  Qi  as  a  parameter.  Solution  for  the  material 

2 

constants  requires  four  equations,  which  are  obtained  by  measuring  T 
for  at  least  four  distinct  values  of  Q^f  where  T2  is  a  real  quantity 
representing  the  magnitude  squared  of  T. 


[t  (f.f'.M'.P"!  V]  2  -  An  n=1’2-3’4 


(98) 


38 


where  ^  is  a  constant  equal  to  the  magnitude  squared  of  the  ratio  of 
the  amplitude  of  the  field  transmitted  through  a  sample  at  angle  Q  to 
the  amplitude  of  the  field  with  no  sample  present.  The  numerical  method 
determines  the  roots  to  the  following  system  of  equations: 

Fn  =  0n)]  2  -  An  =  0  n  =  l,2,3,4  (99) 

The  basic  operation  of  the  computer  code  is  to  calculate  the  values  of 
T  n  corresponding  to  a  trial  vector  of  material  constants,  compare  these 

values  to  the  measured  A^s,  and  revise  the  trial  vector  such  that  the 

2 

values  of  T  *  move  closer  to  the  A  's.  This  process  continues  until  the 
n  n  r 

Fn*s  are  within  a  specified  tolerance  of  zero.  The  criterion  used  to 

2 

compare  the  current  T  n*s  to  the  A^s  is  one  of  least  squares.  A 
function  G  is  defined  as  the  sum  of  squares  of  the  Fn's. 

4  2 

G  =  I  Fn  (100) 

n=  1 

The  code  attempts  to  find  values  of  6  and jj.  that  minimize  G.  It  is 
desired  that  the  numerical  technique  converge  to  a  solution  regardless 
of  the  initial  values.  The  technique  that  was  chosen  is  called  the 
method  of  steepest  descent  (10:511-517).  In  this  method  the  trial 
vector  is  revised  by  evaluating  the  gradient  of  G  at  the  current  values 
of  6  and  fji  and  adjusting  the  trial  vector  in  the  direction  of  tne  nega¬ 
tive  gradient.  The  program  continues  to  revise  the  trial  vector  until: 
(i)  the  value  of  G  is  less  than  the  specified  tolerance 


(ii)  there  is  no  improvement  from  one  revision  to  the  next  (G  does  not 
decrease ) 

(iii)  the  magnitude  of  the  gradient  of  G  is  close  to  zero,  or 

(iv)  the  preset  maximum  number  of  iterations  is  exceeded. 

There  are  other  techniques  that  feature  faster  convergence  than  the 
method  of  steepest  descent,  but  they  require  close  starting  approxima¬ 
tions  to  the  solution.  Although  a  fairly  close  approximation  to  £  an d  jji 
can  usually  be  made  from  measurements  at  other  frequencies,  it  was 
desired  that  the  numerical  solution  not  depend  on  prior  knowledge  of  6 
and  .  If  close  values  are  known  for  the  parameters,  they  can  be  used 
as  the  starting  point  for  the  steepest  descent  technique  to  speed  conve¬ 
rgence.  If  faster  convergence  is  required,  the  steepest  descent  techni¬ 
que  can  be  used  to  obtain  a  suitable  starting  approximation  for  one  of 
the  faster  techniques. 

A  program  utilizing  the  steepest  descent  technique  was  written  for 
both  the  case  of  perpendicular  (MUEPSPERP)  and  parallel  (MUEPSPARL) 
polarization.  In  Chapter  V  of  this  report  the  two  programs  are  tested 

and  the  differences  in  the  results  obtained  are  shown.  The  expressions 

2 

for  T  and  the  gradient  of  G  are  derived  in  Appendix  B. 

Although  MUEPSPERP  and  MUEPSPARL  are  written  for  materials  of 
arbitrary  £  and  jji  ,  it  will  be  shown  in  Section  V  that  it  is  advanta¬ 
geous  to  use  specialized  programs  for  materials  that  are  considered  to 
be  pure  dielectrics  or  for  materials  that  are  considered  to  be  nonmag¬ 
netic.  The  programs  NONMAGPERP  and  NONMAGPARL  are  specialized  versions 
of  MUEPSPERP  and  MUEPSPARL.  The  substitution  jji  has  been  made  and 

T^  is  considered  to  be  a  function  of£?,£  ",  and  the  parameter^  only. 


These  programs  require  only  two  angles  of  incidence  since  there  are  only 
two  unknowns. 

For  the  case  of  pure  dielectrics,  the  only  unknown  is  £  1  and  the 
method  of  steepest  descent  is  not  needed.  The  programs  PURDIELPERP  and 
PURDIELPARL  calculate  a  value  of  £  1  for  each  pair  of  values  Q  and  A. 
These  programs  utilize  a  subprogram  ZEROIN  (11:161-166)  which  determines 
a  root  to  a  function  of  one  variable  by  taking  an  interval  that  contains 
a  root  and  successively  decreasing  the  size  of  the  interval  until  it  is 
smaller  than  a  preset  tolerance  value.  At  this  point  one  of  the  end¬ 
points  of  the  interval  is  output  as  the  root  of  the  equation.  The 
program  must  be  given  a  starting  interval  which  contains  the  root  of  the 
equation.  It  will  be  seen  in  Chapter  V  that  the  equation  for  T  as  a 
function  of  £'  has  multiple  roots.  The  correct  solution  in  this  case 
requires  that  an  interval  containing  only  the  desired  root  be  given  as 
the  starting  interval.  In  Chapter  V  a  method  of  determining  an  appro¬ 
priate  starting  interval  is  discussed. 


V.  Testing  and  Evaluation 


The  measurement  system  is  divided  into  two  parts,  the  computer  code 
that  calculates  £  and  from  measured  values  of  T  and  the  equipment 
that  is  used  to  measure  T  .  Since  comparable  data  for  absorber  mate¬ 
rials  at  94  GHz  is  not  available,  the  two  parts  of  the  system  are  first 
tested  separately  under  known  conditions.  The  codes  are  tested  by 
entering  calculated  values  of  T  for  a  hypothetical  material  of  given 
6  and jJi  and  comparing  the  material  constants  derived  by  the  code  to  the 
given  constants.  The  measurement  of  T  is  evaluated  by  measuring  a 
known  sample  of  plexiglass.  Firally,  results  obtained  for  a  sample  of 
nonmagnetic  absorber  and  a  sample  of  magnetic  absorber  are  presented 
together  with  data  obtained  at  microwave  frequencies  by  other 
techniques . 

Testing  of  Codes 

General  Materials.  The  computer  code  underwent  several  modifica¬ 
tions  and  results  obtained  by  the  different  versions  are  compared.  For 
convenience,  the  following  vectors  are  defined: 

X„:  the  actual  material  constants, 

—a  ’ 

2La  - 

X^:  che  initial  values  chosen  for  the  material  constants, 

Xj-:  the  final  vector  produced  by  the  code, 

A:  the  vector  of  angles  used  in  a  particular  example, 


.  ,  T  •  ’  <4  I  I  M  {  2  ‘  '  ’ 

u  -  ^  ~  _  -  * 

The  inputs  to  the  programs  are  the  test  frequency,  sample  thickness,  A, 
and  T^.  A  separate  code  was  written  for  the  cases  of  perpendicular  and 
parallel  polarization  of  the  incident  wave.  For  the  purpose  of  illus¬ 
tration,  the  results  are  shown  for  a  material  with  constants 
X.  *  [5,  1,  2,  1]  and  thickness  of  100  mils  at  angles 

Q 

A  *  [0°,  20°,  40°,  60°]  and  a  frequency  of  94  GHz.  The  initial  values 

are  chosen  to  be  X*  *  (4,  0.5,  1,  0.5]  .  The  results  are  judged  by 

comparing  the  material  constants  derived  by  the  code  (X^)  to  the  given 

2 

values  (5^)  and  by  comparing  the  values  of  T  (T^  and  T^)  that  corre¬ 
spond  to  and  X^ 

The  results  obtained  by  the  original  codes  are  shown  in  Table  II. 

2 

Although  the  values  of  T  agree  fairly  well  (T^  is  close  to  1^),  the 


derived  material  constants  are  net  close  to  the  actual  ones.  Given 


TABLE  II 


Results  Obtained  with  Original  Codes 


Polarization 

Xf 

If 

T_ 

— m 

G 

4.178 

47.956 

47.960 

0.535 

48.330 

48.318 

8.7  x  10"15 

Perpendicular 

1.927 

49.466 

49.426 

1.297 

51.668 

51.660 

5.041 

47.957 

47.960 

0-569 

48.118 

48.108 

1.0  x  10"14 

Parallel 

1.693 

48.560 

48.524 

_ 

1.120 

49.385 

49.356 

enough  iterations,  the  codes  should  converge  to  X*  but  there  are  practi- 
cal  difficulties  associated  with  using  the  method  of  steepest  descent  to 
solve  Eq  (99). 

Since  T  2S  1  ,  the  error  term  G  can  never  be  greater  than  4  and 

will  typically  be  much  less  than  1  for  lorsy  samples.  For  the  example 
considered  here,  each  element  of  T  is  on  the  order  of  10~-\  When  T  is 

subtracted  from  and  the  difference  squared,  G  becomes  small  quickly. 

2 

To  increase  the  magnitude  of  G  it  was  decided  to  express  T  in  decibel 
form: 


dB 


10  log10(  T  l) 
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TABLE  III 


Results 

Obtained  with 

Final  Codes 

Polarization 

h 

if 

Tm  G 

— m 

4.215 

47.956 

47.960 

0.534 

48.329 

48.318 

Perpendicular 

1.933 

49.460 

8.6  x  10"4 

49.426 

1.294 

51.659 

51.660 

5.064 

47.957 

47.960 

0.547 

48.107 

48.108 

Parallel 

1.839 

48.526 

1.6  x  10-5 

48.524 

1 _ 

1.163 

49.353 

49.356 

Another  difficulty  with  solving  Eq  (99)  is  the  number  of  arithmetic 

calculations  required  to  determine  the  gradient  of  G  (Appendix  B).  The 

number  of  iterations  that  can  be  used  before  G  stops  decreasing  is 

limited  by  the  accumulation  of  round-off  error-  To  reduce  the  round-off 

error  the  codes  were  converted  to  double  precision  FORTRAN.  The  final 

version  of  the  codes  (MUEPSPERP  and  MUEPSPARL)  are  the  double  precision 
o 

versions  with  T  expressed  in  decibels.  The  results  obtained  by  these 
versions  for  the  hypothetical  material  are  given  in  Table  III. 

The  results  have  improved  but  are  still  not  satisfactory.  It  was 
observed  that  MUEPSPERP  consistently  produced  better  solutions  for  £M 
and^jf*  and  MUEPSPARL  consistently  produced  better  solutions  for  £  1  and 
To  demonstrate  this  explicitly,  results  are  given  for  an  example 
with  £  -  fJL  in  Table  IV.  The  reason  for  this  behavior  is  the  different 


1.581 

0.237 

2.084 

0.145 


2.084 

0.144 

1.580 

0.237 


TABLE  V 

Results  Obtained  with  Sequence  of  MUEPSPERP  and  MUEPSPARL 


Perpendicular 

Parallel 

If 

if 

Tm 

if 

T 

— m 

5.066 

47.975 

47.960 

47.975 

47.960 

0.541 

48.333 

48.318 

48.121 

48.108 

1.918 

49.441 

49.426 

48.533 

48.524 

1.185 

51.676 

51.660 

49.365 

49.356 

functional  dependence  of  T  on  each  of  the  four  material  parameters 
for  the  two  different  polarizations.  The  parameter  by  which  the  partial 
derivative  of  T  Is  greatest  dominates  the  gradient  of  G  and  therefore 
the  code  concentrates  primarily  on  improving  the  value  of  that  parameter. 

In  an  attempt  to  improve  the  derived  material  constants  Xf  by  using 


the  equations  for  both  the  perpendicular  and  parallel  cases,  the  corre- 
psonding  codes  were  run  iteratively  with  the  output  produced  by  one 
code  used  as  the  input  for  the  other.  This  approach,  although  cum¬ 
bersome,  yields  a  satisfactory  if-  Also,  if  the  2f  obtained  is  such 

2 

that  the  corresponding  values  of  T  (T^)  closely  approximate  the  meas¬ 
ured  values  (TL,)  for  both  polarizations,  the  confidence  that  the  derived 


constants  are  close  to  X  is  increased.  The  results  obtained  by  running 


the  two  codes  in  sequence  once  are  given  in  Table  V.  X^  can  be  improved 
further  by  continuing  the  iterative  substitution  process. 


Nonmagnetic  Materials.  The  programs  for  nonmagnetic  materials  are 
the  final  versions  of  the  general  programs  with  the  substitution  made 


that  fx  ~jJi0  and  the  derivatives  with  respect  to jj^  and^"  removed. 

n2 


Since  the  measurement  of  T  is  limited  to  four  significant  figures,  it 
is  important  to  examine  the  performance  of  the  codes  when  T^  is  known  to 
four  significant  figures.  The  codes  are  also  tested  by  assuming  a  5% 
error  in  the  measured  values.  The  results  obtained  for  these  cases  are 
summarized  in  Table  VI. 

Pure  Dielectrics.  Although  the  main  application  of  this  system  is 
not  to  measure  pure  dielectrics,  it  was  done  to  verify  the  ability  to 
accurately  measure  the  transmission  coefficient.  This  cannot  be  veri¬ 
fied  with  absorber  materials  directly  because  comparative  data  from 
other  systems  is  not  available  above  26  GHz.  Since  the  codes  developed 
to  determine  £  and^t  for  lossy  materials  do  not  work  well  for  lossless 

n  erials,  a  separate  code  for  pure  dielectrics  was  written. 

2 

T  is  a  function  of  only  one  variable  (£*)  in  this  case,  thus  it 
can  be  graphed  versus  £ 1  for  a  given  frequency,  sample  thickness,  angle 


4  7 


0 

Figure  11.  T""  versus  £’  for  Several  Angles  of  Incidence 

of  incidence,  and  polarization.  Such  a  graph  is  given  in  Figure  10  for 
_/  2 

fc  ranging  from  1  to  10.  It  is  f  »en  that  T  is  not  a  singular  function, 

2 

that  is  , given  a  value  of  T  ,  there  is  not  a  unique  If  one  knows  an 

2 

interval  that  contains  a  unique  value  of  T  ,  then  the  corresponding  £ 1 

2 

can  be  uniquely  determined.  For  example,  in  Figure  10,  T  ■  0.98  has 
seven  roots  for  the  interval  (1,10).  Given  the  following  intervals, 


(1,2),  (2,2.7),  (2.7,4),  (4,4.6),  (4.6,6),  (6,7.1),  and  (7.1,10); 
PURDIELPERP  identifies  the  roots  as  1.543,  2.589,  2.829,  4.529,  4.728, 


6.700,  and  7.185.  A  suitable  interval  can  usually  be  identified  from 
known  results  at  other  frequencies  since  £'  is  not  a  strong  function  of 
frequency  for  lossless  materials.  If  the  measured  T  happens  to  fall 
near  a  critical  point  of  the  graph,  there  will  be  two  possible  roots 
close  to  each  other.  To  determine  the  correct  root,  data  can  be  taken 
at  another  angle  and  the  value  of  £ f  that  is  a  root  for  both  angles  is 
chosen  to  be  correct.  Figure  11  shows  how  the  curve  of  T  versus  £  ' 
varies  as  the  angle  of  incidence  is  changed.  The  critical  points  at  one 
angle  appear  at  different  values  of  £'  at  another  angle.  For  example, 
consider  a  sample  of  thickness  119  mils  being  illuminated  by  a  94  GHz, 
perpendicularly  polarized  waveform.  If  the  measured  values  of  T  are 


0.99 

at 

0 

1.00 

at 

20 

0.86 

at 

40 

and 

0.44 

at 

60 

then  *  2.6  is  the  only  root  that  satisfies  more  than  one  of  the 

measurements.  The  effect  of  measurement  erro’:*  on  the  determination  of  £f 
depends  on  the  instantaneous  slope  of  the  curve  of  T4  versus  £f. 
Therefore,  the  error  in  £*  for  a  given  error  in  T  differs  considerably 
for  different  operating  conditions,  such  as  angle  of  incidence,  sample 
thickness,  polarization,  and  frequency. 


2 

Measurement  of  T 

2 

The  ability  to  measure  T  and  the  influence  of  the  sample-to~ 
antenna  distance  R  were  investigated  by  taking  data  on  a  sheet  of 


TABLE  VII 


Permittivity  of  Plexiglass 


Perpendicular 

Parallel 

ft  R* 

3 

4 

5 

6 

7 

9 

11 

ii 

0 

2.58 

2.60 

2.64 

2.64 

2.61 

2.63 

•..63 

2.60 

5 

2.59 

2.64 

2.69 

2.69 

2.64 

2.62 

2.65 

2.59 

10 

2.61 

2.66 

2.68 

2.71 

2.65 

2.64 

2.67 

2.61 

15 

2.64 

2.69 

2.71 

2.74 

2.68 

2.67 

2.69 

2.64 

20 

2.55 

2.50 

2.47 

2.44 

2.50 

2.52 

2.51 

2.69 

25 

2.60 

2.57 

2.50 

2.52 

2.57 

2.60 

2.58 

2.59 

30 

2.71 

2.65 

2.55 

2.63 

2.65 

2.67 

2.64 

2.64 

35 

2.79 

2.73 

2.59 

2.75 

2.73 

2.74 

2.70 

2.70 

40 

2.92 

2.82 

2.66 

2.87 

2.80 

2.81 

2.74 

2.75 

plexiglass.  The  permittivity  of  plexiglass  is  expected  to  be  very  close 
to  2.6  at  millimeter  wavelengths  (12:334).  Data  was  taken  in  5°  incre¬ 
ments  for  distances  of  3  to  11  inches  and  is  presented  in  Tables  VII  and 
VIII.  The  mean  of  the  data  is  approximately  equal  as  R  is  varied,  but 
the  standard  deviation  decreases  as  R  is  increased  to  11  inches.  It  is 
believed  that  this  can  be  attributed  to  the  greater  interference  between 
the  sample  and  antennae  for  smaller  R  and  to  the  improvement  in  the 
planar  quality  of  the  incident  wave  as  R  increases.  The  sponsoring 
laboratory  has  taken  data  on  plexiglass  on  both  its  time  domain  system 
and  its  network  analyzer  system.  The  data  from  the  two  systems  agree 
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TABLE  VIII 


Hecii  and  Standard  Deviation  of  Plexiglass  Permittivity  Data 


Perpendicular 

Parallel 

R(in)  t 

3 

4 

5 

6 

7 

9 

11 

11 

Mean 

0°-30° 

2.61 

2.62 

2.61 

2.62 

2.61 

2.62 

2.62 

2.62 

Mean 

0°-35° 

2.63 

2.63 

2.60 

2.64 

2.63 

2.64 

2.63 

2.63 

Mean 

0°-40° 

2.67 

2.65 

2.61 

2.67 

2.65 

2.66 

2.65 

2.65 

Std. 

Dev. 

0.34 

0.27 

0.25 

0.36 

0.25 

0.24 

0.19 

0.16 

o 

o 

1 

o 

o 

well;  the  network  analyzer  data  is  presented  in  Figure  12. 

Absorber  Data 

It  was  decided  to  set  R  equal  to  11  inches  for  the  measurement  of 
absorber  and  to  take  data  in  3°  increments  ranging  f :om  9°  to  36° 

The  angles  near  normal  incidence  were  not  used  because  of  reflections 
between  the  sample  and  antennae.  This  choice  yields  10  data  points, 
which  is  sufficient  for  the  numerical  solution.  The  use  of  a  smaller  R 
and  greater  range  of  angles  is  rejected  because  the  reflection  from 
samples  of  absorber  is  typically  greater  than  that  of  the  plexiglass 
sample  and  the  interference  will  be  worsened. 

Nonmagnetic  Absorber.  The  material  chosen  to  measure  was  Eccosorb 
LS-24.  The  sponsoring  laboratory  has  measured  this  material  also  on  the 
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TABLE  IX 
LS-24  Data 


Parallel  Polarization 


Angles 

t 

6 

ff 

€ 

C 

9°, 

U> 

o 

o 

1 

0.951 

0.378 

2.0 

X 

10' 

15°, 

36° 

1.06 

0.401 

1.7 

X 

10' 

12°, 

24° 

0.928 

0.373 

3.5 

X 

10' 

9°, 

24° 

0.929 

0.374 

6.7 

X 

10' 

12°, 

27° 

0.943 

0.376 

4.1 

X 

10' 

15°, 

o 

O 

PO 

0.982 

0.386 

8.2 

X 

10' 

18°, 

33° 

1.02 

0.396 

2.4 

X 

u' 

21°, 

36° 

1.11 

0.415 

2.7 

X 

10' 

18°, 

U> 

o 

o 

1.02 

0.394 

1.3 

X 

10' 

21°, 

33° 

1.04 

0.400 

5.1 

X 

10' 

Perpendicular  Polarization 


9°, 

U> 

o 

a 

1.44 

0.533 

3.4 

X 

10' 

9°, 

21° 

1.49 

0.541 

3.3 

X 

10 

12°, 

24° 

1.51 

0.545 

1.5 

X 

10' 

15°, 

27° 

1.31 

0.510 

1.3 

X 

10' 

18°, 

UJ 

O 

0 

1.56 

0.555 

1.3 

X 

10' 

12°, 

27° 

1.36 

0.518 

2.5 

X 

10' 

15°, 

21° 

1.51 

0.544 

2.5 

X 

10' 

18°, 

24° 

1.53 

0.548 

2  0 

X 

10' 

^  A.  V.  A  AV,  ’  .  -f.  rCm  Ai  *CJ  •<-  A  VVYj  V*  *C.  0-  /. 


X.  The  discrepancy  in  the  results  between  the  perpendicular  and  paral- 

2 

lei  cases  is  believed  to  be  due  primarily  to  measurement  error  in  T 
for  the  perpendicular  data.  The  parallel  data  was  measured  with  a  new 
voltmeter  that  was  more  accurate  than  the  original  meter  used  for  the 
perpendicular  data.  Time  did  not  permit  reconfiguring  the  equipment  to 
remeasure  the  perpendicular  data.  Data  from  the  network  analyzer  system 
is  presented  in  Figure  13. 

Magnetic  Absorber.  The  material  chosen  to  measure  was  FGM-40.  The 
values  of  and  were  arrived  at  by  running  the  pair  of  codes  MUEPSPARL 
and  MUEPSPERP  25  times.  The  final  values  produced  by  the  codes  were 


£  =  19.887  -  j  2.555 

(Parallel) 

fjL  -  0.851  -  j  0.0423 

6  -  19.900  -  j  2.558 

(Perpendicular) 

fj.  =  0.833  -  j  0.00138 

2 

The  values  of  T  (T^)  corresponding  to  the  derived  material  constants 
are  compared  to  the  measured  values  (Tm)  in  Table  XI.  The  codes  took  so 
long  to  converge  because  Jj,  turned  out  to  be  very  close  to jJiQ*  Data  from 
the  network  analyzer  system  was  not  available  for  FGM-40;  data  obtained 
by  the  time  domain  system  is  presented  in  Figure  14  for  comparison  to 


the  results  obtained  here. 


[Qn|i] 


Conclusions  and  Recommendations 


This  thesis  included  the  design  of  both  a  technique  of  measuring 
permittivity  and  permeability  and  the  equipment  needed  to  implement  the 
technique.  A  test  system  was  produced  and  used  to  measure  several 
representative  materials.  This  chapter  presents  conclusions  about  the 
measurement  system  and  makes  recommendations  regarding  further  develop¬ 
ment  of  the  system. 


As  a  result  of  the  experimentation  performed  in  this  thesis,  the 
following  is  concluded: 

1.  The  permittivity  and  permeability  of  lossy,  linear,  homoge¬ 
neous,  isotropic  materials  can  be  determined  from  measurement  of  the 
magnitude  of  the  transmission  coefficient  of  a  planar  sample  at  a  mini¬ 
mum  of  four  angles  of  incidence. 

2.  Data  taken  under  both  perpendicular  and  parallel  polarization 
of  the  incident  wave  is  necessary  to  insure  good  results  for  materials 
of  general  £  and fd. 

3.  Two  data  points  from  either  polarization  are  sufficient  for 
materials  that  can  be  assumed  to  be  nonmagnetic;  the  measurement  system 
performs  well  for  such  materials. 

4.  The  codes  do  not  yield  good  results  for  very  low  loss  mate¬ 
rials,  but  it  is  believed  that  with  further  modification  the  results  for 
such  materials  can  be  improved.  A  separate  code  to  determine  the  per¬ 
mittivity  of  pure  dielectrics  has  been  developed  here  and  yields  satisfac 


tory  results. 


[«] 


5.  The  transmission  coefficient  can  be  measured  at  distances  much 


shorter  than  that  prescribed  by  the  far  field  criterion  by  using  anten¬ 
nae  fitted  with  lenses.  Also;  the  need  to  enclose  the  measurement 
apparatus  in  an  anechoi c.  chamber  can  be  eliminated  by  the  use  of  narrow 
beamwidth  antennae. 

As  a  continuation  of  this  research  effort,  the  following  is 
recommended : 

1.  Since  the  number  of  calculations  needed  for  the  gradient  of  G 
appears  to  be  limiting  the  numerical  solution,  replace  the  exact  expres¬ 
sions  for  the  partial  derivatives  with  a  numerical  differentiation 
routine . 

2.  The  expression  for  the  magnitude  of  the  transmission  coef¬ 
ficient  is  considerably  more  complicated  than  that  of  the  transmission 
coefficient,  yet  the  latter  is  a  complex  quantity.  Since  the  equipment 
has  the  capability  to  measure  the  phase  of  the  transmission  coefficient, 
the  system  could  likely  be  improved  by  developing  a  numerical  solution 
that  solves  for  complex  roots  and  uses  the  phase  measurements. 

3.  As  opposed  to  the  present  approach  of  a  system  of  four  equa¬ 
tions  in  four  unknowns,  nonlinear  regression  may  be  used.  Such  a  tech¬ 
nique  would  have  the  advantage  of  utilizing  all  the  data  points  at  once. 

4.  Other  numerical  techniques  that  feature  faster  convergence  than 


the  method  of  steepest  descent  may  be  investigated.  As  a  first  step,  it 
may  be  desirable  to  use  the  steepest  descent  method  with  one  set  of  data 
points  to  obtain  a  suitable  starting  approximation  and  then  use  a  faster 
technique  for  the  remainder  of  the  data. 
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5.  Since  the  system  now  requires  both  perpendicular  and  parallel 
data  for  general  materials,  it  would  be  an  improvement  to  have  a  single 
code  that  utilizes  both  types  of  data.  One  possibility  is  a  code  with 
four  equations,  two  perpendicular  and  two  parallel,  taken  from  the 
present  codes.  Another  possibility  is  to  use  the  two  present  codes  as 
subroutines  in  a  program  that  would  alternate  between  the  two  subrou¬ 
tines,  using  the  output  of  one  as  the  input  to  the  other.  If  this 
approach  is  followed,  it  would  be  useful  to  examine  the  influence  of  the 
number  of  iterations  within  each  subroutine  in  relation  to  the  number  of 
times  the  main  program  calls  the  subroutines. 

Finally,  modifications  to  the  existing  equipment  could  enhance  the 
performance  of  the  measurement  system.  Specifically,  the  following 
recommendations  are  made: 

1.  Obtain  a  sample  holder  that  can  be  rotated  about  a  horizontal 
axis  as  well  as  a  vertical  axis.  This  would  allow  data  to  be  taken 
under  both  polarizations  without  rotating  the  antennae.  An  alternative 
to  this  is  to  obtain  equipment  that  would  allow  for  the  polarization  to 
be  changed  electrically. 

2.  Develop  the  equipment  so  that  measurements  can  be  made  at  more 
than  one  frequency.  With  the  addition  of  a  network  analyzer  and  a 
computer-controllable  positioner,  the  system  could  be  automated.  Data 
would  be  taken  and  stored  for  each  frequency  at  one  angle,  then  the 
sample  would  be  rotated  and  the  process  continued. 

3.  Experiment  with  increasing  the  antenna  to  sample  distance  to 
beyond  11  inches. 

4.  Fasten  the  sample  holder/positioner  and  waveguide  components 
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rigidly  to  a  support  to  eliminate  vibrations. 

5.  Construct  the  sample  holder  out  of  a  sturdier  material. 

6.  Use  the  original  receiver  equipment  intended  for  this  thesis. 

In  the  course  of  developing  this  material  measurement  system, 
several  unique  characteristics  were  experimentally  verified.  First  of 
all,  measurements  were  made  at  a  frequency  substantially  higher  than 
possible  by  currently  available  systems.  Secondly,  it  was  shown  that 
transmission  could  be  measured  in  free  space  at  very  short  ranges  with¬ 
out  a  shielded  enclosure  by  using  horn  antennae  fitted  with  dielectric 
lenses.  Finally,  a  computer  code  was  written  that  features  global  con¬ 
vergence  to  the  material  constants  when  given  the  magnitude  of  the 
transmission  coefficient  at  four  angles  of  incidence. 

This  thesis  project  has  demonstrated  the  viability  of  the  above 
measurement  approach  and  demonstrated  that  the  approach  is  worthy  of 
further  study.  The  system  can  be  improved  and  extended  in  many  ways, 
however  the  most  pressing  needs  are  to  continue  to  modify  the  computer 
code,  increase  the  range  of  frequencies  at  which  measurements  can  be 
made,  and  automate  the  system. 
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Appendix  A:  Derivation  of  Wave  Impedance  and  Electrical  Thickness 


Figure  15*  Geometry  for  Perpendicular  Polarization 


The  derivations  presented  here  follow  material  given  in  (2:87-90) 


Wave  Impedance 

Perpendicular  Polarization.  The  incident  wave  is  given  by 
x  Ex  ,  y  Hy  ;  where 


Si  =  x  Eix 


Mi  =  y  Hi 
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IX 


e  ^  o 


(103) 


f^)*5  e-(Jko  z) 


(104) 


It  is  convenient  to  define  a  coordinate  system  with  ore  axis  parallel  to 
the  face  of  the  sheet  (r)  and  one  axis  normal  to  the  sheet  (q)-  The 
incident  wave  can  be  considered  to  be  a  uniform  plane  wave  in  the  xy 
plane  propagating  in  the  z  direction  or  as  a  nonuniform  wave  in  the  qr 
frame  propagating  in  the  q  direction-  The  xyz  and  xqr  coordinate  frames 
are  related  by  the  coordinate  transformation 


x  =  x 


y  =  r  cos( ei)  +  q  sin( 6i) 


(105) 


(106) 


=  -r  sin(Ci^)  +  q  cos(0^) 


(107) 


In  the  xqr  frame,  the  incident  field  is  given  by 


=  x  exp 


-  j  *0  C-r  sin(ei)  +  q  cos(0i)]l  (108) 


can  be  determined  from  curl  (E)  =  -jwp.  H 


H.  =  —  curl  (F- ) 

-1  wuo  -1 


(109) 
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curl  (Ei)  -  r  [jj—  (L^)]  “  3  ^Eix^ 


(110) 


since  s  Er  s  0 

Si  *  sfe  -  ^tElx>3]  <111J 

=  ST"  d>j  Ko  cos(ei^  -  ^Cj  *G  sin(e.)]  E.X(112) 

o  L 

=  — -  [r  cos( 6.  )  +  q  sin(?.)]  E.  (113) 

Ci 'll  0  l  1  ix 

Wave  Impedance  of  Free  Space.  Denote  the  wave  impedance 
outside  the  sample  in  the  direction  of  increasing  q  by  ZQ.  Then 


(114) 


(115) 


(116) 


(117) 


sheet,  km,  is  given  by 


8 


The  incident  wave  is  refracted  upon  passage  through  the  air/sample 


boundary  and  propagates  at  an  angle  0T  to  the  surface  normal.  The  form 


of  the  wave  inside  the  sample  is  the  same  as  that  of  the  incident  wave 


with  0^  replaced  by  0r  and  k q  replaced  by  km.  Thus 


+  A  + 

E  =  x  E 
— m  mx 


(121) 


exp  |^-  jk^E-  r  sin(Gr)  +  q  cos(6r)] 


(122) 


H  =  -^  [  r  cos( 9  )  +  q  sin(6  )1  E+ 
— m  (t'u  u  r  ^  r  -  rm 


k  a 

m  _  / ' 


(123) 


A  +  A  + 

=  r  H  +  q  H 

mr  M  mq 


(124) 


where  the  superscript  M  +  "  indicates  the  wave  traveling  in  the  positive 


z  direction  within  the  material  sample.  The  boundary  conditions  require 


that  the  tangential  field  components  be  continuous  for  all  values  of  r 


along  o  ■  0  . 


x  •  Ei  (q  =  0-)  =  X  •  (q  =  0+ ) 


(125) 


(126) 


exp[ j  KQ  r  sin(9i)]  =  exp[ j  km  r  sin(0r)] 


kQ  sin(e^)  =  km  sintG^) 


(127) 


sin(V  =  Nn 

sin(6  )  k 
r  o 


(128) 


([I  c  ) 
Kr  r 


(129) 


Denote  the  wave  impedance  of  the  sample  in  the  direction  of  increasing  q 
by  Zm-  Then 


Z 


m 


,+ 

jnx 

+ 

i 

mr 


(130) 


E 


+ 

mx 


cos(0  ) 
r 


(131) 


05 


m 

- r; - 

(u? )  2  cos ( ) 


(132) 


Z  =  (£)**  sec ( 6  ) 
ir.  t  r 


(133) 


From  Eq  (129) 


sin2( e . )  =  4  €  sin2( e  ) 

i  r  r  r 


(134) 


=  4r  €f  [l  -  cos2  (0r)]  (135) 

68 


(136) 


U  €  -  sin  G.  j 

cos(  e  )  =  [  - i  T2 

“r  €r 


<£)*  ( 4 r  €r)** 

=  [ur  €r  -  sin2(e.)]^ 


(137) 


(tV 

(go); 


Cur  er  -  sin^(e.)] 


(138) 


Normalizing 


Z  /Z 
m  o 


Z  = 


/  u°,  H 

a..  (—) 

x  _ 

2  T£ 
.  €  -  s  in  (  c  •  )  2 

r  r  l  - 


cos ( 6 . 


(^2)1"2 


(139) 


Z  = 


Ur  005(6^) 


2  '  t  o*s 


(140) 


[  u  €  -  sin  (6  •  )n' 

L  r  r  a  “ 

Parallel  Polarization.  Rather  than  proceed  through  a  derivation 
similar  to  that  for  perpendicular  polarization,  advantage  is  taken  of 
the  duality  of  electromagnetic  fields  and  the  result  for  parallel  polar¬ 
ization  is  obtained  by  means  of  the  following  transformation: 


E  ’  =  <£)*  H 


(141) 


H*  =  -  (-)**  E 

—  u 

where  the  primed  variables  denote  the  parallel  case  (2:29). 


(142) 
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(150) 


(£) 


(rr  cos(°r) 


h 


=  (|)^  cos(9r) 


(151) 


[4  €  sin2 ( 0 . )1^ 

LPr  r  1  - 


(a  €  )  2 

r  r 


(152) 


r 


)  2 


€ 


r 


sin2(  9^)~  ‘ 


(153) 


Z  ' 

m 


[a  €  -  sin2(  0  .  )T2 

jr  £  1  J 


€ 

r 


(  154) 


Z  * 


(155) 


(^)**  [^r€r  _  sin2(0i)]'5 

o _ 

40  1, 

€  {—)  *  cos( ei ) 


(  156) 
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ITrtWrt* 


* 


Lur€r  ”  sin“(  ei)-1  ‘ 


€  cos ( c . ) 
r  1 


(157) 


Electrical  Thickness 


6  =  2rr  t 


(158) 


where  T is  the  number  of  wavelengths  inside  the  sheet  along  the  q  direc¬ 


tion  . 


(159) 


6  =  2rr  (■£-) 

Aq 


(160) 


whereA  is  the  wavelength  inside  the  sheet  along  the  q  direction  and  d 


is  the  depth  of  the  material. 


X  .32 

q  kq 


(161) 


where  is  the  wavenumber  corresponding  to\^.  The  propagation  factor 


of  the  wave  inside  the  sheet  along  the  q  direction  is 


»  * 

exp  -  j  K  [-  r  sin(0  )  +  q  cos (6  )] 
m  r  r 


exp(-  j  Krr)  exp(-  j  k^q)  (162) 


Therefore 


k 


q 


k  cos (6  ) 
m  r 


(163) 


:o  (ur€r^  ^r€r  "  sin2(ei^ 

(\i  6  )'2 
r 


(164) 


=  k 


[u  € 
r  r 


.  £  . 
-  sin  ( 


v1- 


(165) 


*  (I2’  ^r-'r 
o 


-  sin2( 


(166) 


6 


(167) 


&  =  (~-^)  C^r^r  -  sin2(ei)]15  (168) 

o 

where\Q  is  Che  free  space  wavelength. 


7  3 


^vi*Qv- *:  « -v:  k  ■£& 


Appendix  B:  Derivation  of  T2  and  the  Gradient  of  G 


Derivation  of 


(Z+l)2  ej5  -  (Z-l)2  e_:i5 


(78) 


_ 4  Z  e~-^ _ 

(Z+l)2  -  (Z-l)2  e"^25 


(169) 


The  superscript  *  is  used  to  denote  the  complex  conjugate  of  a 
quantity . 


T 


★ 

T  T 


where 


* 

T 


A  Z*  eji* 


(Z%1)2  -  (Z*-l ) 2  eJ26* 


(170) 


(171) 


(4  Z  e~ j*  )  (4  Z*  e^5  ) 


[(Z  +  l)2  -  (Z-l)2  e~-'2"  ][  (Z*  +  l ) 2  -  (Z*-l)2  e^26*] 


(172) 


T  2  =  16  Z  Z*  ej(6*“6)  /  C(z+1)2(z*+D2  +  (Z-1)2(Z*-1)2 


J2(6*-5)  -  (Z+1)2(Z*-1 )2  ej2f*  -  (Z*+l)2  (Z-l )2  e-^25]  (173) 
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Since  Z  and  are  complex,  they  can  be  written  in  the  form 


Z  =  a  +  j  b 


6=  s  +  j  t 


1  is  derived  in  terms  of  a,  b,  s,  and  t  and  then  the  expressions 
b,  s,  and  t  are  derived  in  terms  of  £  and 


5*  -  6 .  6*  -  6  =  (s  -  jt)  -  (s  +  jt) 


=  -  j2t 


2  Z  .  Z  Z  =  (a  +  jb)(a  -  jb) 


=  a2  +  b2 


(Z+l)2(Z*f l)2 


VJ 


(Z+1)2(Z*+1 )2  =  [(z+l)(z*+l)]2 

=  [z  Z*+(Z+Z*)  +  i]2 


(  17' 

(17! 

for  a 

(176 

(17- 

( 17{ 

(176 

( 181 

(  181 


Z  +  Z*  =  (a+jb)  +  (a-jb)  (182) 

=  2  a  (183) 

(Z+1)2(Z*  +  1)2  =  (a2+b2+2a -1 )2  (184) 

=  a4+b4+2a2b2+4a( a2+b2+l )  +  6a2+2b2+l  (185) 


(z-i  y 


j 2 ( 5 *-6  ) 


(z-i)2(z*-i)2  eJ2(6*"5)  =  [(z-i)(z*-i)]2  eJ2(6*"5)(186) 
=  [z  Z*  -  ( Z-i-Z  * )  +  l]2  e^2(_J2t)  (187) 

=  ( a2+b2-2a+l ) 2  e4t  (188) 

=  ( a4+b4-*'2a2b2-4a ( a2+b2+ 1  )+6a2+2b2+l  )e4t  (189) 


(Z+1)2(Z*-1)2  =  [(Z+1)(Z*-1)]2 

n  it  ★  ^9 

=  [z  Z  -  (Z-Z  )-l] 


( 190) 

(191) 
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Z  -  Z*  =  (a+jb)  -  (a-jb) 

(192) 

=  j  2  b 

(193) 

(Z+1)2(Z*-1)2  =  (a2+b2-l- j2b)2 

(194) 

(a2+b2-l)2  -  4b2  -  j4b(a2+b2-l) 

(195) 

4 , ,  4 , 2,2  _  2  -.2  ,  •>  ..2  ....  2. .2 

a  +b  +2a  b  -2a  -2b  +1)  -  4b  -j4b(a  +b 

-1 ) ( 196  ) 

4  4  9  2  2  2  22 

+b  +2a  b  -2a  -6b  +l-j4b(a  +b  -1) 

(197) 

*  2  2 
(Z  +1)  (Z-l)  . 

(z*+l)2(z-l)2  =  [(Z+1)2(Z*-1)2]* 

(198) 

4422°2  22 

a  +b^+2a  b  -2a  -6b  +l+j4b(a  +b  -1) 

( 199) 

(Z*+l)2  (Z-l)2  =  m+jn 

(200) 

(Z+l)2  (Z*-l)2  =  m-jn 

(201) 

2  *  2  -i  2  h  * 

( Z+ 1 ;  ( Z  -1)Z  eJZ" 


+  (Z*+1)2(Z-1)2  e~i2° 


=  (m-jn)  e-'2w  +  (m+jn)  e  ^2' 

=  (m-jn)  e-'2^S  +  (nr  jn)  e 

e2*1  [m(e-'2s+e  ^2s)  +  jn(e  -'2s-e-'2s)] 


2t  r  p -^2s  +  e- -^2s 

?  [  2m  (- - - : 


eJ2s  -  e~ -^2s 

2n  - 2j - ^ 


=  2e2t  [  m  cos(2s)  +  n  sin(2s)1 
The  denominator  (Den)  of  Eq  (173)  can  be  written 

Den  =  ( a^+b^+ 2a^b^+4a ( a^+b^* 1 ) +6a^+2b^+ 1 ) 

+  e4t(a4+b4+2a2b2-4a(a2+b2+l)+6a2+2b2+ 
-  2e2t  [m  cos(2s)  +  n  sin(2s)j 


(202) 


(203) 


(204) 


(205) 


(206) 


(207) 


Substituting  for  m  and  n  and  regrouping  terms, 


Den  =  (l+e4t) (a4+b4+2a2b2+6a2+2b2+l)+(l-e4t) (4a) (a2+b2+l) 

-  2e2t[a4+b4+2a2b2-2a2-6b2+ 1 )cos ( 2s ) 

+  4b(a2+b2-l )  sin(2s)]  (208) 

T  2  =  16  (a2*b2Le!!  (209) 

Den 


Derivation  of  a  and  b. 


Perpendicular  Polarization. 


7  _  U  COS  e 

^  ~  “ 7  r 

[u€  -  sin  G"1  ' 


(140) 


where  the  subscripts  "r"  and  i"  have  been  deleted  for  simplicity. 


(tl'-ju")  cos0 


z  = 


9  r 

[  (l. '  -  jji")  (  €  j€")  -  sin  el”5 


(210) 


(u  jn”)  cos© 


[(u'€,-u"€"  -  sin2e)  -  j(€  ' )] 


< 


T-  (211) 
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(220) 


f2  +  g2  =  h[i  +  (a2+&2)"2  -  a  +  (a2+e2)^] 


Therefore 


(a2+e2)^ 


(221) 


c3  =  (a2  +  32) 


cos  e 


[(ji’f+^"g)  +  j(n’g-4"f)]  (222) 


a  =  £22-1  [c ,u'-c2u": 

V2c,  1 


(223) 


b  =  £2=_i  [cjU’-CjH"] 


72  c- 


(224) 


Parallel  Polarization.  From  Eq  (157) 


2  n-lH 


€  -  sin  P 
€  cos  0 


(225) 


cos  0  €’-j€" 


(226) 
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Z  =  - ~ -  [(€'f  +  €"g)  +  j(€"f-€'g)]  (227) 

[(€')+(€")  ]cos  e 


Let 


c4  =  ( € ' ) 2  +  ( €" ) 2 


Therefore 


-3 -  [c1€'+c2€’t]  (228) 

■J2  cos  6 


- - - -  [Cl6"-c2€*1  (229) 

J2  cos  0 


Derivation  of  s  and  t.  from  Eqs  (168)  and  (214) 


6  ■  (7^)  [h  €  -  sin2  el*5 
*o 


(158) 


, 2~d  s  ,  .  N 

U — )  (f  -  jg ; 


(230) 


Therefore 


=  -fl  IT  (A)  Cj 


(231) 


t  =  -  V2  ^(t2-)  c? 

o 


(232) 
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Derivation  of  the  Gradient  of  G 


From  Eq  (100) 


G (  €  '  ,  € " , 4  '  , 4" ) 


4  2 

E  [F  4', 4")] 

n=l 


C 


To  enable  the  use  of  vector  notation,  denote  the  set  of  four  unknown 
material  parameters  ,  and  jj”  by  the  vector 

x  *  [x^,  x2,  x^,  x^]  .  The  gradient  of  G(x),  grad  G(x),  is  given  by 


9radG<*>  = 


(234) 


Denote  any  one  of  the  xn  by  x.  Then 


Thus  tne  derivative  of  Tz  with  respect  to  each  of  the  four  material 
parameters  must  be  obtained.  The  development  will  be  presented  in  terms 
of  a  generic  material  parameter  *'  x  “  as  far  as  possible  and  then 
substitutions  for  the  actual  parameters  will  be  given. 


9  m  2 

Sx  T 


2  _  16  (a^+b^) 

A  +  4B  -  2C  -  8D 


(240) 


4  4  2  2  2  2 

c5  B  a  +  b  +  2a  b  +  6a  +  2b  +  1 


c6  =  (a2  +  b2  +  1) 


-  4  ^  v.4  o  2,2  _  2  ,.2  ,  . 

c-j  =  a  +  b  +  2a  b  -  2a  -  6b  +  1 


cg  =  (a2  +  b2  -  1) 


A  =  c5  (1  +  e  ) 


(241) 


B  =  Cg  a(l  -  e  ) 


(242: 


C  =  Cj  ez  cos (2s) 


(243) 


D  =  Cg  b  e  sin(2s) 


(244 
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and  a  and  b  are  given  by  Eqs  (223),  (224),  and  (228),  (229)  for  perpen¬ 
dicular  and  parallel  polarization  respectively.  Let  "Num"  and  "Den” 
stand  for  the  numerator  and  denominator  of  Eq  (240).  Then 


a  r  m  2n 


Den  (Num)  -  Num  (Den) 
(Den)2 


(245) 


where 


^7  (Num)  =  32e?t  [a^-(a)  +  b^-(b)  +  (a2+fcZ)^(t.)]  (246) 


h  (Den)  *  £(A)  +  4^<B)  -  2^(c)  -  8h(d> 


(247) 


=  4 


4t  a  , .  v 
C5  e  S^U) 


+  ( l+e4t)[a(a2+b2  +  3)^(a)  +  c 


6  bH{fc)] 


(248) 


-4  c6  ae4t  ^(t)  +  (l-e4t)[  (3a2+b2+l)-j^(a)  +  2ab^(b)Z 


(249) 


£<c)-2e2t 


c7[cos(2s)  ^(t)  -  sin ( 2s )  £“(t)] 


+  2  cos(2s)  [c8  a  ^(a)  +  b(a2+b2-3)^(b)n 


(250) 


85 


—  (D)  = 


9x 


cj^2c8b[sin(2s)^(t)  + 


cos(2s)^(s)] 


+  sin  ( 2s )  [2ab^  ( a)+  ( a2+3b2-l  )j^(  b)  ] 


(251) 


The  expressions  for  (30x(s),  (30x(t),  00x(a),  and  30  x(b)  are  derived 
for  each  x  (that  is,  £!,  ,  jjJ  ,  and  jj£)  and  then  substituted  into  the 

above  equations  to  arrive  at  grad  G(x).  The  expressions  for  s  and  t  are 
the  same  under  either  polarization  but  those  for  a  and  b  must  be  derived 
for  each  polarization  separately. 

To  simplify  the  following  expressions  let 

tt  d 


C9 


y[2  X 


K.  =  ip '  +  64" 


-  34 •  -  a4" 


k3  =  a€  ’  +  3€» 


k 4  '  -  a€< 


Perpendicular  Polarization 


&  <*>• 


a  ,  , 

— (s) 


C9  r  K1  ■ 

(  \-  )[u*+(Ti) 


'"1  c3 


is)  = 


cq  ^2 

(  7-  )L-4**+(Ti)] 


(252) 


(253) 


x7 

-  ( — 2^^'  “  CjU")  (265) 

c3 
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-  ( — j )  ( c2^  *  ” 


Parallel  Polarization 


f-  (a)  . 
3x 


3e,(a)  (v'2  cos  -)  cl  +  2c1^ti'+^c3^  +2c2^ 


-  ‘ffXc^  *  C2S”> 


j|^(a)  =  U2  c4  cos  el’1  5§jt-u”-(^>>c2-2^ 


-  I^Hc-E'  +  c,C”) 

c4  1  z 


Trr(a)  =  (72  c,  cos  9) 


i  r  i  ^  o  r  ii 

1  + 


^(a)  =  (72  c4  cos  e)-1  |^[-6.+  (!l)]+|I-c€-+(54): 


(271) 


(b)  • 


1c 

Trr(b)  =  (72  c .  cos  e)"1  -  c7-|— [-u  '  +  (rr) ]+fe- l\i '  +  (r^) ] 


•?  2c  L  4  ' c  2c  'c 

zc2  c3  zcl  c3 


(~” )  (c,  €"  -  c7€«) 
4 


(272) 


at^(b>  *  <V2  c4  cos  6)'1  -  [-u"+(^)] 


-  1^><c2€"  -  c2e'> 


(273) 


e7T( b)  =  (72  c4  cos  e)"1!  -  §^[-€'+(^)]+|£-[€'+(^)]  (274) 


^rr(b)  =  (72  c4  cos  e)_1  -  §^-[€"+(^) (^)]  (27b) 

2  3  1  3 


The  expressions  for  the  nonmagnetic  case  are  not  shown  explicitly.  They 
can  be  obtained  from  the  general  expressions  by  substituting  fj,'  m  1  , 
jjL"  ■  0,  and  deleting  the  derivatives  with  respect  to  jj_'  and 
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A  technique  of  measuring  the  permittivity  and  permea¬ 
bility  of  lossy,  homogeneous,  linear,  isotropic  materials  at 
millimeter  wavelengths  was  developed.  The  technique  v/as 
tested  by  measurement  of  representative  materials  at  a  fre¬ 
quency  of  9 4  GHz. 


For  the  type  of  materials  under  consideration,  the 
permittivity  and  permeability  are  complex  in  general.  Thus 
there  are  four  unknown,  to  be  determined  and  a  minimum  of 
four  measurements  are  inquired.  The  approach  taken  was  to 
measure  the  magnitude  of  the  transmission  coefficient  of  a 
planar  material  sample  in  free  space  at  multiple  angles  of 
incidence.  The  system  of  equations  obtained  from  any  four 
of  the  angles  is  solved  numerically  for  the  permittivity  and 
permeability  of  the  sample. 

The  equipment  consisted  primarily  of  a  Gunn  phase- 
locked  millimeter-wave  oscillator,  a  pair  of  conical  horn/ 
lens  antennae,  a  sample  holder/positioner,  and  a  microwave 
receiver.  Harmonic  mixing  was  employed  to  convert  the  test 
waveform  to  the  operating  range  of  the  receiver. 

y  ■ — 

For  the  case  of  nonmagnetic  materials,  best  results 
were  obtained  if  the  numerical  solution  was  specialized  v.o 
incorporate  the  assumption  that  the  permeability  of  the 
material  is  eaual  to  the  permeability  of  free  space.  For 
a  material  of  arbitrary  permeability,  it  was  necessary  to 
take  data  at  both  perpendicular  and  parallel  polarization 
of  the  incident  wave  to  obtain  satisfactory  results. 
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